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ABSTRACT. We have measured the in situ average electric field attenuation length (La) for 
radio-frequency signals broadcast vertically through the Ross Ice Shelf. We chose a location, 
Moore Embayment, south of Minna Bluff, known for its high reflectivity at the ice-sea 
interface. We confirmed specular reflection and used the return pulses to measure the average 
attenuation length from 75—1250 MHz over the round-trip distance of 1155 m. We find (L a ) 
to vary from ~500 m at 75 MHz to ~300 m at 1250 MHz, with an experimental uncertainty of 
55 to 15 m. We discuss the implications for neutrino telescopes that use the radio technique 
and include the Ross Ice Shelf as part of their sensitive volume. 
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INTRODUCTION 

Several areas of investigation in Antarctica benefit from knowl- 
edge of the attenuation length of radio-frequency transmis- 
' sions in glacial ice. These research areas include radar map- 
ping of subglacial terrain and lakes, location and character- 
ization of subglacial streams, look-ahead radar for over-ice 
traverses, and detection of particle interactions in the ice. 

The authors of this report belong to the community of par- 
ticle astrophysicists aiming to identify high-energy neutrinos 
(produced in astrophysical sources) via the radio emission 
accompan y in g their interaction in Antarctic ice shelves and 
sheets (|l7l ; [3 El S (3) ■ Neutrino interactions initiate a cas- 
cade of charged particles moving in a thin pancake approx- 
imately 0.1 m in extent that traverses tens of meters of ice 
before extinction. The charged particles move near the speed 
of light, c, faster than the velocity of light in ice, which is 
reduced by an index of refraction, n, of approximately 1.78. 
The electric field thus forms an electromagnetic shock wave 
analogous to a sonic boom, known as Cherenkov radiation. 
Given the size of the charged particle pancake, the emission 
among all charges sums coherently for frequencies up to sev- 
eral GHz. This coherence, and the possibility to find high- 
energy neutrino-induced cascades in large natural media such 
as ice, was first identified by Q). 

More recent numerical simulations confirmed and expanded 
the modeling of this emission (| 1 5c 1 2t |25T) . The process was sub- 
sequently confirmed experimentally in a variety of dielectric 
media, including ice, using showers initiated by high-energy 
electron and photon beams dl3 ). 

Because the neutrino detection rate is a strong function of 
the radio-transparency of glacial ice, knowing the attenuation 
length to a reliable precision at several locales in Antarctica is 
of key importance. Extractions of attenuation length for a va- 
riety of Antarctic locations has been ongoing for over several 
decades. For a recent compilation, see Fig. 1 in a previous 
article Q). Most of these values are based on the AC con- 
ductivity of ice (either directly measured or estimated from 
impurity measurements), from which the loss tangent 



and attenuation length are extracted using temperature pro- 
files ifiol : I2H: lioh . In Eqn.[T] e' and e" are the real and imagi- 
nary parts of the complex dielectric constant of the material, 
e = e — ie . From the loss tangent, and its expected vari- 
ation with temperature, the intensity loss in dB m , at a 
frequency u, is 8.686(27r^/2c)(v / e' tan 8). This loss, known as 
the attenuation coefficient, a, is related to the electric field 
attenuation length in meters by 



L a = 1/ln Vl0 a / 10 . 



(2) 



tan 8 — e /e 



(1) 



For a full description of the electromagnetic phenomena of 
radio frequency attenuation in glacial ice see (§) and (0). 

Neutrino searches and the other fields discussed should 
benefit from in situ measurements of the total attenuation 
length, to complement those of the AC conductivity. Several 
such recent measurements have been performed on the ice 
sheet of the Antarctic plateau: one (250-1000 MHz) at Taylor 
Dome in East Antarctica near the Transantarctic Mountains 
(@), one (220-700 MHz) near the South Pole (0), one (140- 
160 MHz) at WAIS Divide (lib , and one using much lower 
frequencies (~5 MHz) at Siple Dome in West Antarctica |27l). 

In this article we report the measurement of the attenua- 
tion length in a floating ice sheet over the range 75-1250 MHz. 
Specifically, we performed a radar bounce through the Ross 
Ice Shelf at a location in Moore Embayment just south of 
Minna Bluff (78°45.022' S, 164°59.291' E). This region was 
chosen because in previous work on basal ice shelf reflectivity 
it was found to be highly reflective to radar (fjjl ). consistent 
with a smooth mirror-like surface below. Having a low-loss 
reflective surface at the bottom allows us to take a measure- 
ment of the returned power and directly extract an attenu- 
ation length. The smoothness of the lower surface of an ice 
shelf (Fimbul Ice Shelf) is confirmed by recent sonar mea- 
surements from the sea below (|24h to be mirror-like down to 
wavelengths of 0.0075 m. The wavelengths considered here, 
within the n = 1.78 bulk ice, range from 0.13 m to 2.25 m, so 
we expect a high reflectivity. For this measurement we assume 
a nominal reflection coefficient of R = 1.0. A true lower value 
would only lead to longer attenuation lengths from our data, 
and for our purpose, R — 1.0 is conservative for neutrino 
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detector design. We show later the effect of reasonable vari- 
ations on R and we find that the effect is small due to the 
presence of the square root and logarithm in Eqn. [2] 

Balloon-borne neutrino observations can include downward- 
pointing neutrino events whose radio emission is reflected 
back upward. It should be noted that for a balloon observa- 
tion of the entire ice shelf, not all locations will have excellent 
reflectivity (|23l) . However, a fixed detector may be placed on 
the surface above the smoothest ice-sea interface. The thick- 
ness of Moore Embayment and its proximity to McMurdo 
station hold promise for this location as a site for a future 
neutrino detector deployed on its surface. Although this ice 
is warmer and expected to have a shorter attenuation length 
than measured on the plateau, our simulations show that it 
will be sufficient to work as a neutrino detector, with an added 
boost from reflected events (Q). The purpose of this work is 
to confirm that the attenuation length is at least as long as 
the value used in that study. Our results indicate attenua- 
tion lengths 1.5-2.5 times longer (better) than used in that 
estimate. 

EXPERIMENTAL APPROACH 

Our basic approach is similar to that of previous work (0;@)- 
We transmit a high voltage impulse, with broadband fre- 
quency content, into the ice and measure the return power 
vs. frequency in another antenna. We compare that to the 
transmission of the same pulse through a short distance, r, 
in air. After accounting for the 1/r 2 geometrical factor, the 
remaining loss is ascribed to attenuation, in which we in- 
clude scattering. We assume 100% reflection at the interface 
of ice and sea water, which yields a conservative value for the 
attenuation length. Systematic uncertainties in the transmis- 
sion and reception are reduced by using only voltage ratios. 
Moreover, because we are measuring over a long distance, i. e. 
several times the attenuation length, the remaining system- 
atic uncertainties in the ratio enter only through a logarithm, 
as shown below. 

The transmitted power per unit area is reduced by two 
effects: a 1/r 2 factor due to the spreading of the wave, and an 
exponential factor exp (— 2r/(L a )) due to losses, where {L a ) 
is averaged over depth. The factor of 2 accounts for power 
being the square of the electric field. We define P v to be the 
power spectral density at frequency v of the pulse as received 
by a 50 receiver and V v to be \fPv x 50S7. We further define 
Vv ice to be V v measured after being transmitted through 
the ice and back to our receiver, and V va i T to be Vv after 
being transmitted horizontally through the air between the 
two antennas separated by d a ; r , typically 9 m. (For the in- 
air calibration, the antenna centers were 2 m above the top 
of the firn. Given the directivity of the antennas, the effect 
of multipath on the calibration is small.) We show later that 
di ce (twice the ice depth) is 1155 m. Hence, before considering 
different transmission efficiency into air versus ice, 

^i/,ice/^i/,air = (^air/^icc) e ^ ^ • (3) 

As described below, we measured T ra tio, the ratio of power the 
antenna transmitted into ice over that transmitted into air, 
accounting for an expected small difference due to impedance 
matching. For electric field measurements, we use %/^ratio- 
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Fig. 1. Unfiltercd output of the high- voltage (HYPS) pulser into 
50 Q. 



Applying this correction twice, once for transmission and once 
for reception, Eqn. becomes 

K,icc/^,air = (Vr^) 2 X (d ail /d icc )e- d ^ /{L -\ (4) 

where T ra ti is a value between 0.8 and 1.2. 

Solving for (La) gives the equation we use to extract the 
attenuation length: 

(L a )=d icc /lnfr ratio ^ ai ^ ah V (5) 

\ ,icc ^icc / 

Note that changes in the transmission into ice versus air, as 
well as any difference in beam pattern, enter only through a 
logarithm. 

Care should be taken not to confuse average attenuation 
length (L a ), which is averaged over depth, with L a . Some 
authors report L a for a particular temperature and depth. 

For our impulse we used a high-voltage pulser (HYPS, 
Grant Applied Physics) designed to drive a 50 load. The 
output of the pulser is shown in Fig. [T] as measured with 
a 3 GHz bandwidth oscilloscope. It output a —2.5 kV peak 
with a 150 ps falltime. We sent and received the impulses 
to the antennas via 23 m Heliax LDF4 50 Q, cables. On the 
receiving end an additional 4 m of LMR240UF cable was 
added. For most data, the signal was amplified by a Miteq 
low-noise amplifier model AFS3-00200120-10-1P-4-L with a 
gain of 49-54 dB from the highest to lowest frequencies. The 
filters discussed below were added after the amplifier and the 
signals were recorded by a Tektronix TDS684 oscilloscope. 

Our main data were taken with quad-ridged horns built and 
designed by Seavey Engineering Associates with a bandpass 
of 200 to 1280 MHz, specifically designed for the Antarctic 
Impulsive Transient Antenna (ANITA) experiment (flil ). Full 
details of the antenna are given in (|13l ). Typical waveforms 
transmitted and received through these antennas are shown 
in Fig. [5] The preferential loss at high frequencies in ice is 
already apparent. 
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Fig. 2. Typical waveforms as transmitted and recorded by the 
quad-ridged horns through 9 m of air (left) and 1155 m of ice 
(right). Both waveforms were recorded without the 900 MHz low- 
pass filter. Each waveform was attenuated or amplified to be ap- 
proximately the same scale on the oscilloscope. 



The transmission band of the antenna can be determined 
from the fraction of power sent to the antenna that is reflected 
back to the pulser, Sn . The Sn for the antenna pointing into 
air (as designed) is shown in Fig. [3] Since the antenna was 
not designed for transmission directly into ice, we made Sn 
measurements in various configurations facing down: directly 
on the snow surface, 1 m above the snow surface, and buried 
0.5 m below the snow surface. We also compared these to 
having the antenna face up to the sky. Little difference was 



Fig. 4. Fraction of incident power reflected back from the Yagi 
antennas (relative to an open termination) while in the snow. The 
two curves are for the two different antennas. 



found among these, and we show Sxi in Fig.[3]for the nominal 
transmission into ice. For lower frequencies we used Yagi an- 
tennas whose Sii are shown in Fig. [4] The critical parameter, 
however, is the ratio of transmitted power: 



T T3 



1 - 10" 



e/10 



1 - 10 Sl 



r/lO' 



(6) 



where Sn is measured in dB. The values of T r atio versus fre- 
quency are shown for the quad-ridged horn and Yagi antennas 
in Fig. [5] Transmission inefficiencies such as Ohmic losses in 
the antenna cancel in the ratios used. 




Fig. 3. Fraction of incident power reflected back from the quad- 
ridged horn (relative to an open termination). The low-pass filter 
began to cut off signal past 900 MHz. 



Fig. 5. T rat i for the quad-ridged horns (black) and Yagi antennas 
(red), with a 900 MHz low-pass filter. Beyond 900 MHz, we use 
1.0 ±0.1. 
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To calculate the distance the radio emissions travel through 
ice, di C e, we measure the two-way transit time between the 
transmitted and received pulse and use a model for index of 
refraction, n, versus depth, z. We measured the total propa- 
gation time, At — 6783 ns, to a precision of 10 ns. We follow 
(|9|) in modeling the ice as two regions: a slab of bulk ice 
with constant n surmounted by a firn layer with varying n(z). 
For consistency with our previous work we take the bulk as 
n — 1.78 ± 0.03, where the uncertainty comes from the range 
of values summarized in (@). 

We use the Schytt model (fir? ) of the firn layer's index of 
refraction, n, versus density, p: 



n(z) = 1.0 + 0.86 p{z), 



(7) 



where p(z) is the specific gravity we measured using core sam- 
ples on the Ross Ice Shelf (Williams Field): 

p(z) = 1.0 - 0.638 e - z / 34 - 7m . ( 8 ) 

In the model, the index of refraction of the firn matches the 
deep ice at z = 67 m, consistent with the results in Fig. 2 of 
(0), beyond which we take n(z) — 1.78. The value of di ce can 
then be found by integrating over the depth: 

o /-djce/2 

At = 6783 ns = - / n(z) dz. (9) 

c Jo 

Even extreme variations in the firn layer model, e.g. linear 
vs. depth, are negligible (<5 m). The dominant uncertainty 
comes from the choice of n for the bulk ice (±0.03), based on 
the distribution of values from (|8|). For our site we measure 
dice = 1155±20 m. This corresponds to an ice depth of 577. 5± 
10 m, agreeing with a measurement of 572 ± 6 m made by 
another team with a different method ~1 km away (ill ). Their 
smaller uncertainty corresponds to a smaller uncertainty used 
on the index of refraction. 

We derive (L a ) from the data using Eqn.[5] The value of V v 
includes the power received in the co- and cross-polarizations 
relative to the transmitted signal. When plotting the data, we 
binned the data points with intervals of 25 MHz, averaging 
Pv,\cc an d -Pi/,air over the width of the bin. The variation 
within 25 MHz was small, but the standard deviation divided 
by the square root of the number of bin elements is included 
in the error bars. To calculate (L a ), we averaged T rat j over 
the width of the bin, although its variation across the bin is 
negligible. 

The final results are shown in Fig. [5] The red data are from 
the Yagi over its frequency band, which was indicated by a 
repeatable and small Sn. The blue data are from the quad- 
ridged horns with a low-pass filter at 900 MHz (Minicircuits 
NLP-1000). The black curve shows a continuation of data 
taken to higher frequencies without the filter up to 1250 MHz. 
We discuss systematic uncertainties and several issues with 
the data below. 

In principle the Yagi data should have also given results up 
to 600 MHz, where we still saw returned power. The T rat ; 
value looked stable to within 20% for this data. However, 
the Six m i° e indicated that up to half the power could 
be reflected and we determined the Yagi data to be unreli- 
able above 200 MHz. Analyzing that data yielded attenuation 
lengths between 200 and 600 MHz about 75 m shorter than 
the nominal quad-ridged horn data. 
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Fig. 6. Calculated average attenuation length, (L a ), as a function 
of frequency. The red line corresponds to data taken with the Yagi 
antennas, and the blue and black lines correspond to quad-ridged 
horn data taken with and without a low-pass 900 MHz filter. The 
modulation vs. frequency is probably an artifact, as described in 
the text. 



At some frequencies we measured a significant amount of 
power arriving in the cross-polarization relative to the trans- 
mitted signal. As shown in Fig. this fraction had a peak 
at about 80% at 450 MHz, which was only seen in ice data, 
not air. This effect could be due either to a rotation of po- 
larization in transit through the ice or due to local resonant 
interaction of the antenna with the ice. In either case, keep- 
ing the received power in the calculation is correct, given the 
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Fig. 7. Fraction of the received signal power that was cross- 
polarized. These data were taken by rotating the receiver by 90° 
after taking the filtered, co-polarized quad-ridged horn data. 
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322. 
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Table 1. Variation of the average extracted attenuation length (in 
m), (L a ), as the assumed power reflection coefficient, R, is varied 
from its nominal value of 1.0. Note that (L a ) only increases for R 
smaller than the nominal value. 



definition of (L a ). We are sure the signal traveled through 
the full depth of ice because of its arrival time. No significant 
delay between the two polarizations was apparent. The posi- 
tion and magnitude of the peak in Fig. [7] did not depend on 
the orientation of the antennas relative to the ice. 

Because of the functional form of Eqn. [5] the main known 
source of systematic error comes from uncertainty in dj ce , 
20 m out of 1155 m, arising mostly from the uncertainty of n 
in the ice. Other sources of error include a 10% uncertainty 
in T ra tio and the variation of power within each 25 MHz bin. 

Our assumption of perfect reflectivity, i. e. reflection coeffi- 
cient R — 1.0, is conservative in that any lower value would 
yield longer attenuation lengths in the analysis. That is, we 
would detect more neutrino events. At our site, the reflec- 
tivity is high (23I) . no worse than R — 0.25 and likely much 
higher. For completeness we show the effect of different as- 
sumed values for 7? in Table [T] Because of the square root 
and logarithm in Eqn. [2] the variation is modest. (|23T I calcu- 
lates a theoretical reflection loss at the seawater/glacial ice 
boundary of —0.77 dB (—0.8 dB if frozen sea ice is present) 
which corresponds to R = 0.83. Had we used that value, the 
effect on our answers would have been longer by typically a 
few percent and never more than 10%. 

An additional uncertainty, probably dominant, is apparent 
in Fig. [6]as a modulation versus frequency of about 183 MHz. 
Tests in the field showed that this modulation depended on 
the local interaction of the antenna with the surrounding di- 
electric snow. The modulation frequency depended on the 
depth at which the antenna was buried or how far above the 
ice it was tested. It also depended on whether the horn was 
filled with snow or just air. However, the modulation could 
not be reduced. Another possibility is that the modulation 
was due to birefringence, by superposing two signals arriving 
at slightly different times due to different propagation speeds 
along the crystal axes. A modulation as observed is consistent 
with a 0.1% difference in wave speeds due to birefringence, 
as has been observed at the South Pole (Q;!!!). Such an ef- 
fect could also be mimicked by partial reflection at or within 
intermediate layers or near the bottom. If due to a single 
layer in the bulk ice, the 183 MHz would correspond to a 
thickness of 0.46 m. Unfortunately, we did not take enough 
data at different angles to confirm or rule out these effects. 
Rather than smoothing out the data, we present the results 



at each frequency bin and assign an uncertainty of approxi- 
mately 55 m at low frequencies and 15 m at high frequencies. 
A three-parameter fit is given below that effectively averages 
out these variations. 

SUMMARY 

We have transmitted and received high- voltage impulses from 
directional horn antennas through the Ross Ice Shelf to the 
sea water below and back. We have demonstrated that the 
reflection at the ice-sea interface is consistent with an ex- 
cellent mirror at these frequencies. In addition, we have ex- 
tracted the average attenuation length over a range of VHF 
and UHF frequencies (75 to 1250 MHz) and found it to range 
from 500 m to 300 m with an uncertainty of order the RMS 
of the variation, or 55 m to 15 m over this range. Since the 
variations are small, if we assume they are an artifact, our 
E-field attenuation length (in m) can be summarized with a 
parameterization over the range 75 to 1250 MHz as 

(L Q ) = 469 - 0.205 v + 4.87 x 10~ 5 v 1 , (10) 

where v is in MHz. The use of a quadratic fit is for convenience 
and not motivated by a particular physical model. The fit val- 
ues, with uncertainties, on the terms in Eqn. [10] are 469 ± 13, 
-0.205±0.036, and (4.87±2.34) x 10~ 5 . Although the data are 
not expected to be described by a constant loss tangent, tan 8, 
the attenuation length decreases with frequency. Naively, the 
data in this range can be summarized by tan 8 ~ 0.0003, 
corresponding to a Q-factor Q = 1/ tan 8 ~ 3300. Few direct 
measurements of the attenuation length have been made over 
these frequencies and our results can be used to constrain 
modeling based on impurities and temperature. 

Our measurements were systematically limited by our un- 
derstanding of birefringence and interaction of the antennas 
with the ice. Based on our experience, we recommend future 
measurements should take data at numerous angles with re- 
spect to the ice fabric. For example, see data taken at Up- 
stream B Camp (|19l ). results which warn that the crystal 
fabric can change orientation on scales of 100 m. Care should 
also be taken in antenna design or simulation to understand 
the interaction with the dielectric medium. The impulses were 
strong enough that reflected signals from intermediate layers 
were visible, so a mapping of birefringence vs. depth should 
be possible if many angles are taken. 

These attenuation lengths are long enough, comparable to 
the depth, that neutrino detectors proposed for this region of 
the Ross Ice Shelf appear viable. In addition, balloon flights 
that have observed this region are also sensitive to neutrino 
interactions. The coherence of the pulse is well maintained 
with no distortion beyond that due to the attenuation. The 
lengths are conservative in that if the mirror at the ice-sea 
interface is not perfect, the actual attenuation lengths would 
be longer. 
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